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We calculate within a mean-field theory the spectral signatures of various striped d-wave super- 
conducting phases. We consider both in-phase and anti-phase modulations of the superconducting 
order across a stripe boundary, and the effects of coexisting inhomogeneous orders, including spin 
stripes, charge stripes, and modulated density- wave. We find that the anti-phase modulated 
(i-wave superconductor exhibits zero-energy spectral weight, primarily along extended arcs in mo- 
mentum space. Concomitantly, a Fermi surface appears and typically includes both open segments 
and closed pockets. When weak homogeneous superconductivity is also present the Fermi surface 
collapses onto nodal points. Among them are the nodal points of the homogeneous d-wave super- 
conductor, but others typically exist at positions which depend on the details of the modulation 
and the band structure. Upon increasing the amplitude of the constant component these additional 
points move towards the edges of the reduced Brillouin zone where they eventually disappear. The 
above signatures are also manifested in the density of states of the clean, and the disordered system. 
While the presence of coexisting orders changes some details of the spectral function, we find that 
the evolution of the Fermi-surface and the distribution of the low-energy spectral weight are largely 
unaffected by them. 

PACS numbers: 74.81.-g, 74.25.Jb, 74.20.-z, 74.72.-h 



I. INTRODUCTION 



Over the last decade it has become increasingly clear 
that the cuprate high-temperature superconductors ex- 
hibit inhomogeneous electronic structures both in their 
"normal" and superconducting states^. In particu- 
lar, several scanning tunneling spectroscopy experiments 
have produced evidence for real-space inhomogeneity 
of the superconducting gap without apparent spatial 
ordering^^, while others have demonstrated periodic 
modulations of the local density of states (LDOS)^^^. 
Even if much of the modulated signal can be attributed to 
interference patterns resulting from scattering of quasi- 
particles off impurities^' at least part of it is likely 
associated with a spatially periodic structure of the su- 
perconducting order and possibly of other coexisting elec- 
tronic orders^^^^. 

In this paper we study, within mean field theory, var- 
ious realizations of modulated d-wave superconducting 
(DSC) phases. We focus on stripe phases in which the 
lattice-translational invar iance is broken in one direction, 
and calculate their spectral signatures in the clean limit 
and in the presence of disorder. The quasiparticle den- 
sity of states of clean translational-symmetry breaking 
states has been calculated by Podolsky et al. in Ref. 20. 
There, the effects of a small modulated order parameter 
added to a large- amplitude d-wave superconductor, were 
investigated. Here we extend this treatment to strong 
modulations, consider also the distribution of spectral 
weight in momentum space, and study additional sce- 
narios which were not treated in Ref. 20. In particu- 
lar, we analyze the case of an anti-phase modulated su- 
perconductivity, where the (i-wave order parameter suf- 
fers a TT-phase shift across each stripe boundary. Such 
a state has been proposed recently^^'^^ as the source for 



the apparent decoupling between the Cu-0 planes in 1/8 
doped La2-ccBaa;Cu04^^. Similar structures were also 
found in a mean-field study of the DSC resonating va- 
lence bond phase of the t — J model^^. Finally, we supple- 
ment the extensive literature on the Fourier-transformed 
local density of states (FT-LDOS) of disordered super- 
conducting phases, with and without coexisting inhomo- 
geneous orders^^^^, by calculating the FT-LDOS of a 
system with anti-phase modulated superconductivity. 

We demonstrate that in sharp contrast to its ho- 
mogeneous and in-phase modulated counterparts, the 
anti-phase modulated d-wave superconductor exhibits a 
Fermi surface which includes extended parts of the non- 
interacting Fermi-surface, as well as closed pockets. The 
corresponding zero-energy spectral weight appears pre- 
dominantly over arcs in momentum space, whose extent 
shrinks with increasing amplitude of the order parameter. 
We find that this general behavior is robust with respect 
to changes in the details of the band-structure and in the 
functional form of the modulations. Furthermore, the 
general evolution of the Fermi-surface and the distribu- 
tion of low-energy spectral weight is largely unaffected 
by the introduction of additional coexisting orders. The 
situation changes significantly when a homogeneous su- 
perconducting component is also present in the system. 
Under such conditions the Fermi-surface collapses onto 
nodal points, which include those of the homogeneous d- 
wave superconductor, and generically, additional points 
at positions which depend on details of the modulation 
and band-structure. Upon increasing the constant order 
parameter the extra nodal points move towards the edges 
of the reduce Brillouin zone (BZ), where they eventu- 
ally disappear. The presence of the Fermi-surface is also 
reflected in the low-energy density of states of the anti- 
phase modulated DSC state, and in the momentum-space 
structure of the FT-LDOS of this system when impurity 
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scattering is included in the analysis. 



II. THE MODELS 

Our starting point is a tight binding model of electrons 
hopping on a square lattice. In the following we con- 
sider two non-interacting band structures, corresponding 
to Bi2Sr2CaCu208+(5 and La2-a?Baa;Cu04. While most 
of the scanning tunneling experiments have been carried 
out on the first, owing to the good quality of surfaces 
which can be obtained in this system, the latter is prob- 
ably the most promising compound for the observation 
of the signatures we discuss in the following, especially 
in the context of anti-phase modulated superconductiv- 
ity. Consequently, most of our results will be presented 
for this material. To model Bi2Sr2CaCu208+5 we use the 
tight-binding Hamiltonian provided by Norman et. al.^^ ^ 
for the system near optimal doping. The 1/8 hole-doped 
La2-a:Baa;Cu04 is described by the free dispersion^^ 
^(k) = — (ti/2)[cos(/Ca;) + cos(/c^)] — ^2 cos{kx) cos{ky) — jJ.^ 
where ti = 1.72 eV, ^2 = — 0.15ti and /i is the chemical 
potential which we adjust in order to maintain the re- 
quired level of hole-doping. 

In this work we study, using a mean-field Hamiltonian, 
the spectral signatures of a spatially modulated (i-wave 
superconductor. We also consider within the same ap- 
proximation, the effects of other inhomogeneous coex- 
isting orders on such properties. Motivated by the ex- 
perimental evidence for the existence of "stripes" in the 
high-temperature superconductors^, we analyze the case 
where the discrete translation symmetry of the under- 
lying lattice is broken by the periodic orders in the x 
direction. Since often the observed charge and spin mod- 
ulation periods are 4 and 8 lattice constants, respectively, 
we assume a unit cell of either lengths for the modulated 
orders. In the following we describe the real-space struc- 
ture of the various configurations which are studied in 
the paper. 



A. The superconducting order 

We consider two configurations for the DSC order. In 
the first its magnitude is sinusoidally modulated across 
each site-centered stripe, but its phase is constant over 
the sample. This "in-phase" configuration is described 
by the mean-field Hamiltonian 
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where x and y are measured in units of the lattice con- 
stant which we take from now on to be 1. Here, and 
throughout the paper = 7r/4. Note, that since the or- 
der parameter does not change its sign across the stripe 
boundary the resulting configuration has a 4-site unit 
cell, as presented in Fig. 1. On a microscopic level it 
may refiect a variation in the pairing amplitude as a re- 
sult of a corresponding modulation of the hole density in 
the stripe phase, as discussed below. 
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FIG. 1: The DSC configurations considered in this work: a) 
site-centered in-phase modulated, and b) site-centered anti- 
phase modulated, c) bond-centered anti-phase modulated. 
The magnitude of the order parameter is depicted by the 
width of the lines, while its sign is given by the colors of 
the bonds. 

The second, "anti-phase", configuration is similar to 
the first, except that the phase of the order parameter 
changes by tt across a stripe boundary. It has been sug- 
gested in Ref. 21, that despite being unconventional, 
such a negative Josephson coupling between stripes may 
occur in 1/8 doped La2-xBaa^Cu04 , with significant con- 
sequences for the transport through this system. We con- 
sider the site-centered version of the configuration, see 
Fig. 1, and analyze the Hamiltonian 

H^-dSC = ^ < ^cos[g^(x + 1/2)] 
x,y I 



[^x,y]'-x^l 

^ cos{qxx) 



yi '-x,yi'-x^l,y] 



-i,y]^i,y+li ^i,yi^i,y+l] 



H.c. . (2) 



H.c. 



(1) 



The anti-phase nature of the configuration results in an 
8-site unit cell. 

We have also considered the bond-centered analogs of 
the above site-centered configurations. Fig. 1 depicts, for 
example, the bond-centered anti-phase modulated DSC. 
We have found, however, that the results are hardly af- 
fected by this detail, and therefore in the following we 
report them for the site-centered orders. 
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B. Coexisting electronic orders 

Various electronic orders have been proposed to ex- 
ist in the cuprate high-temperature superconductors, es- 
pecially in the underdoped regime. In this paper we 
consider the case where the modulated superconducting 
order is accompanied by one or more of the following: 
charge stripes, anti-phase spin stripes and an anti-phase 
modulated (i-density wave (DDW). 

0©®®0©®©0 
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FIG. 2: A sketch of the magnetic and charge distribution in 
the stripe model. The arrows represent the amplitude and 
sign of the spin density, and the shading of the circles repre- 
sents the charge density (the darker the circle is, the higher 
the charge density). 



For the charge stripe phase we assume a site-centered, 
period-4 sinusoidal charge density {ux^y^ + ^x,yi) = 
<I^CDW |sin(g'a^x)|. When spin stripe order is also present 
the boundaries of the charge unit-cells serve as anti-phase 
domain walls for an antiferromagnetic (AF) spin order of 
the form {rix^y} - rix^yi) = (j)AF{-lT^^ s'm{qxx). The 
combined charge-spin stripe configuration is depicted in 
Fig. 2, and the corresponding Hamiltonian is described 
in appendix A. 

The uniform DDW state was proposed as a model for 
the pseudogap state in the underdoped cuprates^^. It 
describes a condensation of electron-hole pairs with non- 
zero angular momentum into a phase with staggered cur- 
rent loops. In the following we consider the anti-phase, 
period-8, modulated version of this state as displayed in 
Fig. 3. Note that current conservation has been incor- 
porated into the model. For a detailed expressions of 
the order parameter and the mean field Hamiltonian see 
appendix A. 




FIG. 3: A schematic representation of the modulated DDW 
order parameter. The arrows represent the currents in the 
system. Gray arrows carry half the current of the black ones. 



III. RESULTS 

A. The Fermi-surface and the distribution of 
low-energy spectral weight 

It is well known that the Fermi-surface of a homo- 
geneous d-wave superconductor consists of four nodal 
points. A very different behavior emerges in the case 
of an anti-phase modulated DSC, as shown by Fig. 4. 
The figure exhibits the evolution of the Fermi-surface of 
La2-xBaa;Cu04 with increasing strength. A, of the anti- 
phase modulated DSC, where the chemical potential was 
adjusted to maintain a 1/8 hole-doping level for each set 
of parameters. One can clearly see that the model sup- 
ports an extended Fermi-surface. As the amplitude A 
increases, segments of the non-interacting Fermi-surface 
near (0,7r) and symmetry related regions, become pro- 
gressively gapped. At the same time the Fermi-surface 
develops closed pockets which shrink with increasing A 
and eventually disappear completely. The process pro- 
ceeds until a single pocket remains when A is of the order 
of the band-width. It then continues to shrink upon fur- 
ther increase of the DSC order. This behavior is generic 
as can be seen in Fig. 5, where similar results are pre- 
sented for Bi2Sr2CaCu208+5 . The rate at which the 
Fermi-surface evolves, depends, however, on the band- 
structure. The fiatter sections in the anti-nodal regions 
of Bi2Sr2CaCu2 08+(5 tend to develop a gap more quickly 
than the rest of the Fermi-surface. 




FIG. 4: The Fermi Surface (a)-(c), and the low-energy spec- 
tral weight (d)-(f) of a model corresponding to 1/8 hole-doped 
La2-xBaccCu04 in the presence of an anti-phase modulated 
DSC. The Fermi surface is plotted in the first quadrant of the 
Brillouin zone. The spectral weight is integrated within a 20 
meV Lornetzian window centered at zero energy. The ampli- 
tude of the modulated DSC is A = in (a,d), A = 0.075ti in 
(b,e) and A = 0.15ti in (c,f). 

We believe that the origin of the Fermi-surface can be 
traced to the formation of zero-energy Andreev bound 
states on the anti-phase domain walls of the DSC or- 
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at points which are related to them by the modulation 
wave vector. 
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FIG. 5: The Fermi Surface and the low-energy spectral weight 
in a model of 17% hole-doped Bi2Sr2CaCu2 08+(5 with anti- 
phase modulated DSC. The spectral weight was integrated in 
a similar manner to Fig. 4. The results are given for A = 
in (a,d), A = 0.075t in (b,e), and A = 0.15t in (c,f) (|t| 0.6 
eV, see Ref. 38.) 



der. A single junction between two phase-biased d- 
wave superconductors has been studied by Tanaka and 
Kashiwaya^^. They have calculated the quasiparticle 
LDOS at the interface between the superconductors and 
demonstrated that for a 7r-phase shift the LDOS displays 
a pronounced zero-bias peak, associated with zero-energy 
Andreev bound states. The latter move to higher ener- 
gies when the phase difference is tuned away from tt, and 
correspondingly the LDOS exhibits a pseudo-gap behav- 
ior at low biases. In our model we consider a chain of 
such junctions and therefore must take into account the 
multiple scattering processes between them. We have not 
carried out a detailed analytical analysis of this problem. 
Nevertheless, the numerical real-space structure of the 
zero-energy wavefunctions reveal that the states residing 
near the end-points of the open Fermi-surface segments 
are indeed localized on the domain walls. On the other 
hand, the zero-energy states on the Fermi-pockets and 
away from these end points are approximately evenly dis- 
tributed over the system. 

Angle resolved photoemission spectroscopy (ARPES) 
is arguably the most suitable probe for detecting the ex- 
tended Fermi-surface. More generally, ARPES measures 
the spectral function A(k, uj) of the system. In Figs. 4, 5 
we plot the distribution of low-energy spectral weight of 
an anti-phase modulated d-wave superconductor, namely 
the frequency integrated A(k, cj) within a 20 meV Lor- 
netzian window centered at zero energy. As shown in the 
figures, the spectral weight lies predominantly along arcs 
in momentum space. As we increase the superconduct- 
ing order the arcs shrink towards the BZ diagonals. For 
the sinusoidal modulation studied here the distribution 
of spectral weight along the arcs is continuous. However, 
for a square-wave modulation, gaps open on the arcs at 
the boundaries of the reduced BZ at kx = =t7r/8 and 




FIG. 6: The momentum occupation function n(k) of a model 
of 1/8 hole-doped La2-xBaxCu04 with an anti-phase modu- 
lated DSC of strength A = 0.075ti. 
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FIG. 7: The evolution of the nodal points in a model of 1/8 
hole-doped La2-ccBaccCu04 with anti-phase modulated DSC 
of strength A = 0.075ti, as function of the amplitude, Ao, 
of an additional homogeneous DSC component. The curve 
depicts the Fermi-surface for the case Ao = 0. 

ARPES can also provide a measurement of the momen- 
tum occupation function n(k) = J^^{du/27r)A{k,ij). 
It is this function which depicts most directly the pres- 
ence of the closed segments, or "pockets", in the Fermi- 
surface. Fig. 6 displays n(k) for the anti-phase mod- 
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FIG. 8: Fermi surfaces (a)-(d) and low-energy spectral weight distributions (e)-(h) of a model corresponding to 1/8 hole- 
doped La2-ccBaxCu04 with anti-phase modulated DSC and additional orders. The amplitude of the DSC order parameter is 
A = O.lti, and the additional orders are: in (a,e) CDW with (i)CDW — 0.075ti, (b,f) spin-stripe order with (\)af — 0.075ti, 
(c,g) charge and spin stripe order with (\)af — (j^CDW = 0.075ti, (d,h) modulated DDW with (I)ddw = 0.0125ti. Note that the 
reduced BZ is halved in the presence of spin stripes or a modulated DDW, thus complicating a straightforward comparison 
between the different Fermi surfaces. 



ulated DSC in La2-a:Baa^Cu04 . While Fermi-surface 
pockets in the normal state can always be classified as 
electron-like or hole-like, this is no longer true for the 
superconductor whose quasiparticles are linear combina- 
tions of both electrons and holes. However, as can be seen 
from Fig. 6, the pockets dit kx — and kx — ±37r/4 are 
close to being electron-like, while the pockets at kx — 
±7r are of the hole-like type. Recently, much interest 
has been generated by the observation of Shubnikov-de 
Haas and de Haas- van Alphen oscillations in underdoped 
YBa2Cu307_5 in strong magnetic fields which suppress 
superconductivity^^'^^. The experimental results sug- 
gest that the Fermi-surface contains coherent electron- 
like pockets. Theoretical studies of non-superconducting 
spin-charge modulated phases have found evidence for 
electron^^ and hole^^ pockets. Although unlikely to be 
relevant to the understanding of the aforementioned ex- 
periments, it would be interesting to investigate the im- 
pact of the pockets which we find in the anti-phase model, 
on the transport and magnetic properties of the system. 
We, however, did not pursue this direction. 

When a constant DSC component, of arbitrary small 
magnitude, is added to the anti-phase modulated sys- 
tem the Fermi-surface changes dramatically and collapses 
onto nodal points. Among them we always find the nodal 
points of the uniform d-wa^e superconductor along the 
diagonals of the BZ. However, as demonstrated by Fig. 
7, additional nodal points typically appear for not too 
large magnitude, Aq, of the constant component. For 



the parameters used by us their spectral weight is about 
a half of the weight of the points along the diagonals. 
When Aq is increased, the additional points move from 
the tips of the pockets, where they originate, towards 
the edge of the reduced BZ where they disappear, while 
the points along the diagonals stay put. The number 
and positions of the nodal points vary as function of 
the strength of the modulated order and of the band- 
structure, but the general behavior remains the same. In 
the case of an in-phase modulated DSC we always find 
a Fermi-surface which consists of the four nodal points 
of the uniform superconductor. This is not surprising 
since the in-phase modulated order can be represented 
approximately as an equal strength combination of anti- 
phase and constant components, for which, as seen in 
Fig. 7, the extra nodal points have long vanished. It is 
worth mentioning that while the change in the topology 
of the Fermi-surface takes place even for an infinitesimal 
Aq, the original Fermi-surface of the anti-phase modu- 
lated system is only gapped by an amount of order Aq. 
Therefore, the low-energy spectral weight distribution of 
a system with weak additional homogeneous DSC is close 
to the one presented in Figs. 4, 5. 

We have also studied the Fermi-surface and the low- 
energy spectral weight in the presence of additional elec- 
tronic orders, as described in section II. At large, we 
have found that the coexisting orders do not change the 
general features outlined above. The Fermi-surface of 
the in-phase modulated system is still comprised of the 
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conventional wave nodal points. The introduction of 
the orders into the anti-phase modulated superconductor 
does have some effect on the Fermi-surface and the spec- 
tral weight distribution, as seen in Fig. 8 (Note that the 
reduced BZ is halved in the presence of spin stripes or 
a modulated DDW.) Most notably, we observe that the 
CDW order induce gaps on the Fermi-surface at multiples 
of the ordering wave vector. Although the shape of the 
Fermi-surface is somewhat different, its extent is largely 
unaffected by the additional orders. This is also true for 
the general evolution of the Fermi-surface with increasing 
A and with the addition of a uniform DSC component. 




Energy (meV) 



B. The density of states 

The low-energy signatures discussed above are also 
manifested in the LDOS, which can be measured using 
scanning tunneling spectroscopy. Fig. 9 presents the den- 
sity of states of the modulated superconducting states af- 
ter averaging over the position in the sample. The energy 
dependence of the DOS in the case of the in-phase modu- 
lated DSC is close to that of a uniform d-wave supercon- 
ductor whose amplitude is given by the zero-wavevector 
component of the modulated order. In particular, the 
DOS vanishes at zero-bias (In order to simulate the finite 
energy resolution of the experiments we have introduced 
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FIG. 9: Averaged density of states in a model of LBCO with 
(a) anti-phase modulated DSC and (b) in-phase modulated 
DSC, for various magnitudes of the superconducting order. 



FIG. 10: LDOS in a model of La2-xBaccCu04 with an anti- 
phase modulated DSC of strength A = O.lti. Atom is 
located on the symmetry axis of the stripe - the left most site 
in Fig. 1. 
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FIG. 11: P(7r/2,o) in a model of LBCO with (a) anti-phase 
modulated DSC and (b) in-phase modulated DSC, for various 
magnitudes of the superconducting order. 



a 5 meV Lorentzian broadening into our numerical cal- 
culation. This is the reason for the apparent zero-energy 
DOS in the figure.) In contrast, the zero-energy DOS 
of the anti-phase modulated system remains finite. We 
find that the low-energy behavior of the DOS depends on 
the band-structure. For La2-xBaajCu04 the DOS con- 
tains a small dip at zero bias whose width does not scale 
with the strength of the modulated superconducting or- 
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der, see Fig. 9. In Bi2Sr2CaCu208+(5 on the other hand 
(not shown here), there is a more pronounced dip which 
does scale with A. In both cases the low-energy DOS 
decreases when A is increased, in agreement with the re- 
duction in the extent of the Fermi-surface, as described 
above. 

Although the unit-cell in the anti-phase superconduct- 
ing system is composed of 8 sites, the real space period- 
icity of the LDOS is halved, as shown in Fig. 10. This is 
a result of the symmetries of the model. Translating the 
system by 4 lattice constants is equivalent to the change 
A —A, which does not affect the LDOS, as long as the 
average pairing potential is zero. An additional constant 
DSC component breaks this symmetry, opens a d-wave 
gap in the LDOS, and reverts its period to 8 sites. Conse- 
quently, the first non- vanishing Fourier component of the 
LDOS, {pqiuo))^ in both the in-phase and the pure anti- 
phase modulated models is at q = (7r/2,0). Its energy 
dependence is shown in Fig. 11. While for the in-phase 
model it exhibits a symmetric structure around zero bias 
within the low- volt age regime, such symmetry is absent 
in the anti-phase case. 

The addition of electronic orders to the superconduct- 
ing system does not change the main features described 
above. Zero-bias spectral weight is still observed in the 
anti-phase DSC model and is absent in the in-phase sys- 
tem. The periodicity of the LDOS is also unaffected by 
the CDW and the DDW orders, while the spin-stripe or- 
der introduces new Fourier components at multiples of 
q = (7r/4,7r). 



C. FT-LDOS 

In recent years, scanning tunneling spectroscopy mea- 
surements of the FT-LDOS have become a much-used 
tool in the quest for identifying unconventional electronic 
orders in the cuprates^^. We have calculated the effects of 
quasiparticle scattering off a single impurity on the FT- 
LDOS in a system with modulated superconductivity. To 
this end we have implemented the T-matrix approxima- 
tion, as described in appendix B. 

Since most of the spectral changes in the modulated 
system occur at low energies we choose here to concen- 
trate on the zero-energy FT-LDOS, and present, in Fig. 
12, its distribution for an anti-phase modulated DSC sys- 
tem with, and without, an additional homogeneous DSC 
component. As can be seen, the Fourier transformed 
spectra is highly detailed. However, this fine structure is 
unlikely to be accessible experimentally as it is of very 
small weight compared to the prominent peaks which 
dominate the distribution. Consequently, we concentrate 
on the latter. 

In both cases studied by us the strongest peaks appear 
at multiples of the modulation wavevector q = (7r/4,0). 
[Note that disorder eliminates the symmetry which pro- 
duces the period-4 LDOS of the clean anti-phase modu- 
lated system, thus enabling the appearance of the peaks 
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FIG. 12: FT-LDOS at u; = of a model of LBCO with 
(a) anti-phase modulated DSC order and (b) coexisting anti- 
phase and homogeneous DSC orders. The amplitude of the 
modulated order is A = O.lti and that of the constant compo- 
nent is Ao = O.Olti. The results are for an impurity strength 
of Vb = 150 meV. Darker regions correspond to higher FT- 
LDOS and we have indicated high-intensity peaks by circles: 
Peaks associated with the modulation wave vector are encir- 
cled in blue. Peaks that correspond to scattering between the 
nodal points of the homogeneous DSC are marked in green. 
Peaks due to scattering between the additional nodal points 
which appear in the mixed phase are marked in pink. Other 
high- intensity peaks which are not simply related to scatter- 
ing between nodal points are designated in red. 



at (±7r/4,0).] The positions in k-space of many of the 
remaining peaks may be associated with wavevectors 
connecting high zero-energy LDOS regions of the clean 
system^' This is especially true for the case with addi- 
tional homogeneous DSC where the zero-energy spectral 
weight appears at isolated nodal points. We were able 
to identify, as shown in Fig. 12, peaks which are asso- 
ciated with scattering between the nodal points of the 
homogeneous d-wave superconductors, as well as peaks 
which originate from scattering between the extra nodal 
points that are induced by the presence of the modulated 
DSC order, see Fig. 7. It is more difficult to trace the 
root of the peaks in the FT-LDOS of the pure modulated 
system, since it possesses an extended Fermi-surface. 



IV. CONCLUSIONS 

We have identified several signatures which set apart 
the anti-phase modulated (i-wave superconductor from 
other d-wave superconducting states. Most notably, we 
have shown that it supports a continuum of zero-energy 
excitations which form an extended Fermi-surface. The 
existence of such a Fermi-surface requires a pure 7r-phase 
shift of the superconducting order across the domains 
boundaries. The presence of a constant DSC order, no 
matter how weak, drastically modifies the manifold of 
zero-energy states and results in a Fermi surface which 
consists of isolated nodal points. Nevertheless, much of 
the low-energy spectral weight still concentrates along ex- 
tended arcs in momentum space, and can be detected by 
ARPES measurements. ARPES can also be used to iden- 
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tify the additional nodal points which we predict exist 
under such circumstances, beside the conventional nodal 
points associated with a homogeneous d-wdive order. The 
periodic modulation of the superconducting order can be 
traced using scanning tunneling spectroscopy, which can 
also detect the accompanying low-energy spectral weight, 
and the identifying features in the FT-LDOS of the anti- 
phase modulated superconductor. 
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APPENDIX A: THE HAMILTONIANS 



1. Superconducting orders 



The real-space representation of the in-phase modu- 
lated superconducting Hamiltonian is given by Eq. (1). 
Since only integer x and y are of interest, when convert- 
ing it to the momentum-space representation, we approx- 
imate the absolute value function by its first eight Fourier 
components, with the result 



/o(k) = — 2V2 cos{qx/ 2) cos kx — {I -\- V2) cos ky 
8 L 

A r ^ 

/2(k) = -- v2sin(g^/2)(cos/c^ - sin/c^) - cos/c^ 

o L 

A 



/4(k) = -{-1 + V2) COS ky, 

fe{k) = /2(k + 6q). 



(A5) 



In the presence of an anti-phase modulated supercon- 
ducting order one obtains 



/i(k) = — [cos{kx + - cos ky] , 



/7(k) 



[cos{kx - qx/2) - cos ky] . (A6) 



The analysis of the various orders was carried out 
using mean-field Hamiltonians written in momentum 
space. The modulated superconducting orders are 
conveniently expressed in terms of the spinor t/^J^ = 

(4t' • • • ' c-ki, c_(k+q)|, . . with q the ordering 

wavevector, and k varying over the reduced Brillouin 
zone (RBZ) associated with the order. In this basis the 
Hamiltonian 



i^k V^k 



\<ieRBZ 



takes the form 



vAk Ck 



^k 



(Al) 



(A2) 



where .4k is a diagonal matrix with entries ^(k),^(k + 
q), . . ., and where Ck contains the particle-particle cou- 
plings. 

While the RBZ of the in-phase modulated supercon- 
ductor is twice that of its anti-phase counterpart, we 
choose to describe both in the RBZ of the latter, i.e., 
— 5a;/2 < kx < (?a;/2 and — tt < ky < it, where q = 
(7r/4, 0). In this case Ck is an 8 x 8 matrix, 



Ck 



/ /o(k) 
/r(k) 

V/7*(k) 



/i(k) 
/o(k + q) 



/7(k) \ 
/6(k + q) 

/o(k + 7q) J 



(A3) 



characterized by the functions /n(k). 

The only non-vanishing contribution to the uniform 
superconducting order is 



/o(k) 



Ao 



(cos kx — cos ky). 



(A4) 



2. Non-superconducting orders 

For the Hamiltonians of the non-superconducting or- 
ders we use a similar representation to Eq. (Al), but with 
V^k defined with respect to the wavevector q = (7r/4,7r). 
Here, we choose to let k run over the RBZ of the mod- 
ulated DDW and spin-stripe orders, despite it being half 
the size of the RBZ of the charge-stripe phase. The result 
is 



^^k 



^k + Sk 

-A - ^Sk 



(A7) 



where 77 = 1 for the charge and DDW stripe phase, and 
= — 1 in the case of spin-stripe order. Sk contains the 
particle-hole couplings 



Bk 



/ 50 (k) 51 (k) 
91* (k) 5o(k + q) 

V57(k) 



57(k) \ 
96(k + q) 

(;o(k + 7q) J 



(AS) 



The non-vanishing contributions in the charge-stripe 
phase are 

5o{k) = 5£^[l + 2cos(g,)], 



52 (k) 
54 (k) 



4 

56 (k) = 

^CDW 



CDW 



[1 -2cos{qx) 



(A9) 



where the absolute value function was again approxi- 
mated using the low harmonics. The ^o(k) term renor- 
malizes the chemical potential. The latter was chosen to 
maintain the desired hole doping. 
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The spin-stripe phase is characterized by the foUowing 
terms 



^3(k) = -^5(k) = i 



(AlO) 



We considered an anti-phase, period-8, modulated ver- 
sion of the d-density wave, (see Fig. 3). The correspond- 
ing real-space order parameter is 



{ 

(All) 

where (x) is a periodic variation of the usual heavy- 
side function, i.e. 

{1 2Tin/qx < X < (2n + 1)tt / Qx 
-1 (2n - l)Ti/qx <x< 27rn/qx (A12) 
X = utt/qx 

for all integers n. The corresponding functions entering 
the momentum-space Hamiltonian are given by 

^i(k) = /i(k,5g,/2), 

^3(k) = /i(k,7g,/2), 

^5(k) = /i(k,-7g,/2), 

gr{k) = /i(k, -5g,/2) , (A13) 

where 

h{k,qx) = ^Dwe'^"sin(4g^)[l + 2cos(2g^)] 

X [cos(A:^ + Qx) - cos{ky) cos(g^)]. (A14) 

3. Combining the orders 

When superconducting and non-superconducting or- 
ders coexist we need to expand our basis to include them 
both. We use 

'^k = (4T'4+(g,,0)T' • • • '4+(7g,,0)T'4+(gx,7r)T' • • • ' 
^k+(7ga,,7r)T' ^-ki' <^-k+(g^,0)i5 • • • , C_k+(7g^,7r)i)- 

(A15) 



The combined Hamiltonian is 32-dimensional and is de- 
fined in the momentum-space region —qx/2 < kx < qx/2, 
and — 7r/2 < ky < 7r/2. The terms in this Hamiltonian 
are determined by adding the corresponding terms, con- 
necting the same momenta, in Hamiltonians (A2) and 
(A7). 

APPENDIX B: CALCULATING THE FT-LDOS 



The clean system can be described in terms of the 
spinor ipl. defined in appendix A. In this basis the Hamil- 
tonian is H = XlkGi?BZ ^k^kV^k, where is a model 
dependent matrix. The Green's function of the clean 
system, G^(k, cj), is obtained by analytically continuing 
G^(k, iujn) = [i(^nl — ^k]~"^, via iujn ^ OJ -\- where / 
is the identity matrix. In the numerical calculations we 
have used an energy broadening of S =0.5 meV. 

The impurity scattering problem can be simplified by 
working in the Nambu basis in real space, (clpCri). The 
transformation of the Green's function between the two 
representations is given by 



GO(r,r';a;) 



1 



(27r)2 



E 



^ik(r-r') 



(Bi) 



l^eRBZ 

(/>t(r)G?i(k,a;)0(rO 0t(r)G?2(k, u;)0(rO 
0t (r)GOi (k, u;)0(rO 0t (r)GO, (k, u;)0(rO 



where (/)t(r) = (1, e^^^ e^^^^ . . .), and G^-, ij = 1, 2, are 

the four blocks comprising G^{k,uj). 

We have included the impurity scattering in the T- 
matrix approximation^^, which implies the following 
equation for the electronic Green's function 

G{r,r';iOn) = G\r,r';io) (B2) 

+ J dYidY2G\Y,ri-Uj)f{vi,r2\L0)G\v2,v'-L0). 

For simplicity we have assumed scattering off a non- 
magnetic ^-function impurity centered at the origin. Un- 
der such conditions the T-matrix is given by T(r, r'; cj) = 
[Uo"V3 - G^{{)]uj)]-^6{v)6{y'), where Vb is the impurity 
strength and r are the Pauli matrices. The local den- 
sity of states at a point r in the sample is then given by 
/)(r, uj) = — ;^Im[Gii (r, r; uo) + G22 (r, r; —uj)] , from which 
the FT-LDOS is readily obtained. 



^ For a review of some of the experimental evidence for the 
existence of meso-structure in the cuprates and its possible 
relation to the mechanism of high temperature supercon- 
ductivity see E. W. Carlson, V. J. Emery, S. A. Kivel- 
son, and D. Orgad, in "The Physics of Superconductors: 



Superconductivity in Nano structures, High-Tc and Novel 
Superconductors, Organic Superconductors'^ Vol 2, p. 275, 
edited by K. H. Bennemann and J. B. Ketterson (Springer- 
Verlag 2004). 
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